Spin-triplet supercurrent carried by quantum Hall edge states through a Josephson 

junction 
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We show that a spin-polarized Landau level in a two-dimensional electron gas can carry a spin- 
triplet supercurrent between two spin-singlet superconductors. The supercurrent results from the 
interplay of Andreev reflection and Rashba spin-orbit coupling at the normal-superconductor (NS) 
interface. We contrast the current-phase relationship and the Fraunhofer oscillations of the spin- 
triplet and spin-singlet Josephson effect in the lowest Landau level, and find qualitative differences. 

PACS numbers: 73.23.-b, 73.43.-f, 74.45. +c, 74.78.Na 



I. INTRODUCTION 

The coexistence of the quantum Hall effect with the 
superconducting proximity effect provides a unique op- 
portunity to study the flow of supercurrent in chiral edge 
states. The usual quantum Hall edge states^ in a two- 
dimensional (2D) electron gas are created by the inter- 
play of cyclotron motion and reflection from an elec- 
trostatic potential, propagating in a direction dictated 
by the cyclotron frequency lo c — eB/m. At the inter- 
face with a superconductor Andreev reflection from the 
pair potential takes over, converting electrons into holes. 2 
Since the sign of both the effective mass m and charge 
e change upon Andreev reflection, the cyclotron rotation 
keeps the same direction for electrons and holes and the 
chirality of these Andreev edge states is preserved^ - — 

While the superconducting proximity effect is short- 
ranged in the direction perpendicular to the edge states, 
it is long-ranged in the parallel direction. Indeed, a su- 
percurrent can flow through a 2D electron gas even if the 
magnetic field is so strong that only a single Landau level 
is occupied — provided the spin splitting by the Zeeman 
effect is sufficiently small; 7 -! 8 - Andreev reflection from a 
spin-singlet superconductor couples opposite spin bands, 
so spin polarization of the Landau level suppresses the 
supercurrent j2»1£1 

Recent studies of ferromagnetic Josephson junctions 
have shown that a spin-triplet proximity effect (with elec- 
trons and holes from the same spin band) can be induced 
by a spin-singlet superconductor, if the spin is not con- 
served at the ferromagnet-superconductor interface^ - — 
In the 2D electron gas of a quantum well formed in a nar- 
row band gap semiconductor, such as In As or InSb, the 
Rashba effect is a significant source of spin-orbit coupling 
in quantum Hall edge states.— When contacted with Nb 
electrodes, these structures show a strong proximity ef- 
fect in the quantum Hall effect regime.— ~— 

In this article we investigate whether the spin-polarized 
lowest Landau level of a 2D electron gas can carry a 
spin-triplet supercurrent between two spin-singlet super- 
conductors, as a consequence of the Rashba effect on 
Andreev edge states. We find that a long-range spin- 
triplet proximity effect does exist, with a critical current 



oc (d/l so ) 2 , determined by the spin-orbit scattering length 
Z so in the normal region and the distance d over which the 
electrostatic potential drops upon entering the supercon- 
ductor. It is a small effect, but the fact that it exists as 
a matter of principle opens up the possibility to optimize 
it. 

We calculate the current-phase relationship (depen- 
dence of the supercurrent on the superconducting phase 
difference) and the Fraunhofer oscillations (dependence 
on the magnetic flux through the junction) of the spin- 
triplet Josephson effect and compare with the corre- 
sponding spin-singlet effect. Some of our spin-singlet re- 
sults are know n 7 ' 8 ' 18 , but some are new. In particular, 
we find a complete suppression of the Fraunhofer oscil- 
lations in the spin-singlet case for a critical value of the 
width W of the Josephson junction. (These spin-singlet 
results may be of interest also for graphene, which shows 
a strong proximity effect^ without significant spin-orbit 
coupling.) 

In Sec. |TT] we formulate the problem of edge state trans- 
port along a superconductor, in the form of an effec- 
tive Hamiltonian in the lowest spin-split Landau level. 
The parameters entering into this Hamiltonian are de- 
rived from the Bogoliubov-De Gennes equation in the 
Appendix. The spin-triplet Josephson effect is analyzed 
in Sees. IIIII and IIVI and compared with the spin-singlet 
counterpart in Sec. [V] We conclude in Sec. I VII 



II. SPIN-POLARIZED TRANSPORT ALONG A 
SUPERCONDUCTOR 

A. NS interface 

We consider the scattering by a superconductor (ex- 
citation gap A , Fermi energy Ep,s ^> Ao) of a single 
spin-polarized edge channel in a 2D electron gas in a per- 
pendicular magnetic field B. The lowest Landau level at 
\hhj c ± ig/iB-B is split by the Zeeman energy g^B, and 
spin polarization is ensured by taking the Fermi level Ep 
in the 2D gas in between the two spin-split levels (typi- 
cally Ep s» Tjhuj c ). 

The characteristic energy and length scales at the 
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FIG. 1: Schematic drawing of the energy scales and length 
scales at an NS interface. The electrostatic potential profile 
is shown as a blue solid curve, the Fermi level is a red dashed 
line, and the superconducting excitation gap is green dashed. 
The red solid lines indicate the spin-split lowest Landau level, 
with only a single spin band occupied (short black arrows). 
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FIG. 2: Dispersion relation of edge states along an NS inter- 
face in the lowest spin-polarized Landau level, for the electron- 
like mode (solid) and the hole- like mode (dashed). The black 
curves are calculated in the Appendix from the Bogoliubov-De 
Gennes equation (for g^B.B -C huj c , va/v c <C 1, \/l m "C 1). 
The red lines are the small-p approximation ()2.6|) . 



normal-superconductor (NS) interface are shown in Fig. 
[TJ On the superconducting side we have the coher- 
ence length £o = ^f,s/Ao, and the Fermi wave length 
Xp,s — ^/kF.s — wfiVFjS/EFjS- We require that £o is 
small compared to the magnetic length l rn = ^h/eB, to 
ensure that B is well below the upper critical field of the 
superconductor . 

The electrostatic potential step at the NS interface ex- 
tends over a distance d, which we assume to be inter- 
mediate between Xp,s an d Co- These length scales are 
therefore ordered as 

A F:S «d«£ <Z m . (2.1) 

We include the rounding of the electrostatic potential 
step because it has a major effect on Andreev reflection. 
(For an abrupt interface, d < Aj?, Andreev reflection is 
strongly suppressed even without spin polarization.) The 
step in the pair potential is also rounded, but this has 
no significant effect on Andreev reflection (since Ao -C 
Ef,s)- 

On the normal side of the NS interface the Fermi wave 
length is l m , so this is not an independent length scale. 
The spin-orbit scattering length and coupling energy are 
l so = h 2 /ma and _E so = ma /H , respectively, with a the 
Rashba coefficient. 



B. Edge channel Hamiltonian 

The wave function \I/ = (ip e ,iph) of the electron and 
hole excitations (both in the same spin band) is an 
eigenstate of the Bogoliubov-De Gennes Hamiltonian H 
with energy eigenvalue e (measured relative to the Fermi 
level). Electron- hole symmetry dictates that if (ijj e ,iph) 
is an eigenstate at energy e then (ip^, tp*) is an eigenstate 
at energy — e. This requires 



where the Pauli matrix acts on the electron-hole degree 
of freedom. 

At low excitation energies an effective Hamiltonian, 
containing only terms linear in momentum along the 
edge, is sufficient. The form of this effective Hamiltonian 
is fully constrained by the requirements of Hermiticity 
and electron-hole symmetry, 

rr_ 1 f{v c ,p- eA} {v A ,p} \ r9 
2\ {v* A ,p} {v c ,p + eA})- W 

Here s and s are coordinate and unit vector along the 
edge, p = —id/ds is the canonical momentum, and A = 
As is the vector potential in a gauge where it is parallel 
to the edge. (We set H = 1 in intermediate formulas and 
write +e for the electron charge.) The anti-commutator 
{a, b} = ab + ba ensures that H is Hermitian even if the 
velocities v c and i>a depend on s. 

The gauge transformation <F t— > exp(^xcr z )^ , trans- 
forms the Hamiltonian as follows, 

H ^ e ixa *He- ix ' 7 * = 

l({ Vc ,p-eA- X '} {\vA\e^ +2 ^, P }\ . . 
2 \{\v A \e-^- 2 ^,p} { Vc ,p + eA + X '})' 1 ' 

with x' = dx/ds and v A — |«A|e^. We ensure that 
v A is real positive by chosing 2% = —(j>. The effective 
Hamiltonian then takes the form 

H = (v c + v A a x )p - ev c Aa z - \i{v' c + v A a x ). (2.5) 

C. Dispersion relation 

For s-independent A, v c , and v A the momentum p 
along the edge is conserved. The Hamiltonian (|2.5p de- 
scribes two chiral modes with dispersion relation 



cr x H*a x = -H, 



(2.2) 



e = v c p ± \J (ev c A) 2 + (v A p) 2 , 



(2.6) 
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see Fig. [2j At e = the two modes have the same group 
velocity u group = de/dp, given for v& <C v c by 



"A 



/v c . 



(2.7) 



Let us express v c and v\ in terms of the character- 
istic parameters of the NS interface. As derived in the 
Appendix, the two velocities v c and va are given, up to 
numerical coefficients of order unity, by 



v c d 



(2.8) 



The velocity v c is the same as the cyclotron drift velocity 
in the lowest Landau level along a normal, not super- 
conducting boundary, in the limit of a steep confining 
potential. The confinement by the superconductor is ef- 
fectively in that limit because the penetration depth £o 
of the edge state into the superconductor is less than its 
transverse extension l m . 

The velocity f A which governs the coupling of electrons 
and holes is smaller than v c by a factor d/l so . Although 
it is the superconducting order parameter which scatters 
electrons into holes (Andreev reflection), the dependence 
on Ao drops out in the regime £o ^ ^m- The ratio d/l so 
appears in the calculation in the Appendix as the product 
of two factors, with a cancellation of the magnetic length: 
One factor is the probability d/l m of Andreev reflection 
with change of spin band and the other factor is the spin- 
flip probability Z m /Z so - The length Z so refers to spin-orbit 
scattering in N. There may also be spin-orbit scattering 
in S, but that would contribute to «a a much smaller 
amount of order v c (d/l so )(d/l m ) 2 , see the Appendix. 



D. Effect of screening current 

The vector potential along the NS interface is deter- 
mined by the screening of the magnetic field from the in- 
terior of the superconductor^ Consider an interface at 
y = with the superconductor in the region y < 0. The 
edge state propagates in the +x direction. The vector po- 
tential is A = A(y)x, with magnetic field B = —A'(y)z. 
We denote by A = A(0) the value of A at the NS in- 
terface. The Andreev-Rashba edge channel extends over 
a distance l m from the interface, so the effective value of 
the vector potential is Aar — Aq — c m l m B. The value of 
c m s» 0.88 is calculated in the Appendix. 

The value of Ao follows from the London equation for 
the screening supercurrent density j, 



1 



J 



2e/ioA 2 \ds 



2eA Q 



(2.9) 



with A the London penetration depth. For l m > A the 
magnetic field decays exponentially oc e~' y '/ A upon en- 
tering the superconductor. (This is the Meissner phase of 
a type-II superconductor, reached for magnetic fields be- 
low the lower critical field.) The screening current within 



a distance £o <C A from the interface is j = B/hqA. In 
the gauge where the order parameter is real, one thus has 
A = -XB. 

We conclude that the vector potential A in the edge 
state Hamiltonian (12.51) takes the value 



(2.10) 



4ar = —(c m l m + A)B, 



along the NS interface in the Meissner phase l m > A. 
The phase difference 27t|Aar|/Vo accumulated per unit 
length between electron and hole (with ipo = h/2e the su- 
perconducting flux quantum) is increased by the screen- 
ing current. This is a Doppler effect of Andreev reflection 
from a moving superconducting condensate i 20 ' 21 

For magnetic lengths in the range £o < lm < A the 
magnetic field penetrates into the superconductor in the 
form of Abrikosov vortices. In this mixed phase j4ar 
depends on the detailed configuration of vortices. We 
will consider this regime by treating ^4ar as a random 
function of the position along the NS interface. 



E. Transfer matrix 

We transform from a Hamiltonian to a scattering de- 
scription of the edge channel transport, which is the 
description we will use to calculate the Josephson cur- 
rent in a superconductor-normal-metal-superconductor 
(SNS) junction. 

The particle current operator is 



J = dH/ dp — v c + v&<j x 



(2.11) 



We require < «a < v c , so J 1 / 2 is real Hermitian. To 
construct a unitary transfer matrix we transform H to 

H = J- l / 2 HJ- 1 ' 2 =p- J-^eVcAazJ- 1 / 2 . (2.12) 



The wave function '5 = J 1 / 2 ^ then satisfies 



(2.13) 



The transfer matrix M(s2,si) relates the function 
\P(s) at two points along the boundary, \&(s2) = 
M(s2, Si)^(si). Integration of Eq. (|2.13|) gives the ex- 
pression 



M(s 2 ,si) 
= V a exp 

= V s exp 




(eJ- 1 + J- 1/2 ev c Aa z J- 1/2 
e(v c - va<J x ) ev c Aa z 



(2.14) 



with V s the operator that orders the noncommuting ma- 
trices from left to right in order of decreasing s. 

The transfer matrix is unitary, M' — M~ l , as an 
expression of particle current conservation: ( V E'| V 1 / ) = 
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FIG. 3: Left panel: Superconducting ring, enclosing a mag- 
netic flux $, interrupted in one arm by a normal segment 
(non-shaded region, containing a flux 5&). Right panel: en- 
largement of the SNS junction between the normal (N) and su- 
perconducting (S) regions. The normal region is a 2D electron 
gas in the quantum Hall effect regime, with a spin-polarized 
edge channel near the Fermi level (dashed, with arrows indi- 
cating the direction of motion). 



at point s1 ut = {-W/2, 0+) on the left edge. The scat- 
tering matrix for this process is S\(e) = M(s° ut , s^)\ e . 
Similarly, a wave incident on interface 2 from point 
sJf = (— W/2, _ ) on the left edge comes out at point 
s° ut = (W/2,0~) on the right edge, with scattering ma- 
trix S 2 (e) = M(s° ut , 4 n )| e . 

The SNS junction is a segment of a ring enclosing 
a magnetic flux accounted for by a vector potential 
A$> = $S(y)y (for \x\ < W/2). The total scattering ma- 
trix S(e) for a closed scattering sequence, from s 1 / 1 to s° ut 
to s 2 n to S2 Ut to s™, is given by 

S(e) = e"' J ^ /{f ' S2{e)e- na ^ /vo S l {e). (3.1) 

The contribution to the scattering matrix from the nor- 
mal segments of the boundary can be calculated imme- 
diately from Eq. (|2.I4I) . because va — and no operator 
ordering is required. We thus obtain 

S(e) = e ^o e i(Vvo)(*+«5*/2) ( T, M2 ( £ ) 

x e -i{*/<po)(*-s*m*, Ml ( s y (3,2) 



('I' | J I'll/) is independent of s. Electron- hole symmetry 
is expressed by 



Ml 



(2.15) 



Since the expression (|2 . 14[) does not assume that the 
parameters v c ,va are uniform along the edge, it may 
also be used to describe the transport along a boundary 
containing both normal and superconducting segments. 
On the normal segments va = (no electron-hole cou- 
pling), while the cyclotron drift velocity v c is still given 
by Eq. (|2.8I) (for a confining potential that is steep on 
the scale of l m ). The vector potential A along the nor- 
mal edge is determined by the enclosed magnetic flux, 
without the correction (|2.10l) from the screening current 
that is present along the superconducting edge. 

Consider a superconducting segment connecting two 
normal boundaries. An electron enters the supercon- 
ducting segment at the left end and exits at the right 
end, either as an electron or as a hole. At e = the 
transfer matrix M commutes with a z . This implies that, 
at the Fermi level, the electron exits as an electron with 
unit probability. At finite excitation energy e the proba- 
bility for Andreev reflection (from electron to hole, with 
the transfer of a Cooper pair to the superconducting con- 
densate) vanishes as e 2 when e — > 0, in accord with Refs. 



The flux through the junction is <5$ = BLW and tq = 
§ dsv c /(v 2 — v A ) « 2(L + W)/v c is the time it takes a 
quasiparticle to circulate along the entire perimeter of 
the junction. The matrices M n give the contribution to 
the scattering matrix from the interface with S„ (without 
the scalar factor, which has already been accounted for 
in the factor e lST °): 



M n {e) = P s exp 



ds 



-ev A cr x 



ev c Aa z 



. (3 - 3) 

The Josephson current /( ( f ) ) flowing in equilibrium at 
temperature T through the SNS junction is related to the 
scattering matrix by2^ 

1 A 00 

/ ( $ ) = 2 d$ 51 2k * T ln dct t 1 - ( 3 - 4 ) 

p=0 

The imaginary energies are Matsubara frequencies, iuj p — 
(2p + l)iirkBT. The prefactor 1/2 accounts for the fact 
that only a single spin band contributes to the supercur- 
rent. (In Ref. [23J it is canceled by the spin degeneracy.) 
The Josephson current is a periodic function of the flux 
<!> through the ring, with period ipg. The critical current 
I c of the Josephon junction is the largest value reached 
by |/($)|. 



III. EDGE CHANNEL JOSEPHSON EFFECT 

The geometry of the SNS Josephson junction is shown 
in Fig. [3] It consists of two parallel NS interfaces, inter- 
face 1 at y = L/2 and interface 2 at y = —L/2 (for both 
interfaces \x\ < W/2). A wave incident on interface 1 
from point s 1 / 1 = (W/2, + ) on the right edge comes out 



IV. SPIN-TRIPLET SUPERCURRENT 

A. Calculation 

To calculate the supercurrent we use the fact that 
va/vc is a small parameter. An expansion in this pa- 
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rameter is made possible by the identity 



V a exp 



ds [a(s) + b(s)} 



w 



A(W)V s exp[ dsA- 1 (s)b(s)A(s) , (4.1) 



A(s) — V s ' exp ds a 



(*') > 



(4.2) 



valid for any pair of operator functions a(s), 6(s). An 
easy way to prove this identity is to call the right- 
hand-side X{W) and calculate dX/dW = [a(W) + 
b(W)]X(W). Integration then produces the left-hand- 
side. 

With the help of Eq. (|4~Tj) , the expression ([53)1 for the 
scattering matrix M n along the interface with S n takes 
the form 



M n (e) =e M ^^exp 

, ev c A n (s') 



w 

is I ds 
o 



9 2 



t/n(s) 



ds' 



(4.3) 
(4.4) 



The integral in the definition of the phase U n {s) runs over 
a distance s along the NS„ interface, and a n is the total 
phase accumulated by the vector potential A n (s) along 
that interface. 

To first order in va the expression (|4.3p reduces to 



M„(e) = e M "^ - iee M "^ 



a x I ds(v A /v, 
Jo 



2\ e 2iU n a z 



1 -ieAJ 
-ieA„ 1 



with the definitions 



(4.5) 

(4.6) 
(4.7) 



From Eq. f|3.2f) we obtain the determinant, to second 
order in va, 



f V A / f S 

A„ = / ds — exp ( 2i / ds' eA n {s') 
7o w c V Jo 



ds eA n {s). 



Det [1 - S{iw)} = 2e 



cosh(wTo) 



cos(7rJ$/^ + a x + a 2 ) - ie^^O A x | 2 + |A 2 | 2 ) 



; 2 Re (/ Xl / x * e ^2~ ai +2^/ Va ) 



(4.8) 



and then substitution into Eq. (I3.4[) gives the supercur- 




-2-10 1 2 

8$/tp + ip'/ir 

FIG. 4: Low-temperature critical current I c as a function of 
the flux 5$ through the normal region, plotted from Eq. (|4.11[1 
in units of e|AiA2|/ro . 



rent, 



2irk B T 
Wo 



x u? v [cosh(o;pro) - cos(7n5$/<^ + ol\ + a 2 )] . 

(4.9) 

This expression holds for arbitrary temperature and 
for arbitrary variation of A(s), v c (s), and va(s) along 
the two NS interfaces, which is fully accounted for by the 
parameters A„ and a n [Eqs. (|4.<j[) and (|4.7[> ]. We will 
now discuss this general result in some illustrative limits. 



B. High and low-temperature regimes 

The high-temperature limit {UbTtq ^> 1) of Eq. 
is given by the p = term in the sum over Matsubara 
frequencies, 

/($) = 47r 2 esin(27r$/^o + ^)(fc B T) 3 |AiA 2 |e- 7r ' SBTTo . 

(4.10) 

The low-temperature limit is obtained by replacing the 
sum by an integration, with the result 

/($) = - l Al f 2 l sin(27r$/^ + iP)T{ir6<5>/<p + 

IT T Q J 

(4.11) 

ip = axgAi — argA 2 + a 2 - a%, ip' = otx + a 2 , (4.12) 

r°° to 2 

J r (x)= du . (4.13) 

J cosh lu — cos x 

The function F(x) oscillates between T(0) = 2ir 2 /3 and 
J-(tt) = vr 2 /3. 

The vector potential along the NS interfaces introduces 
a phase shift if) in the sinusoidal current-phase relation- 
ship, as a result of which the current /($) is no longer 
and odd function of the flux $ through the ring. 

In the high-temperature regime the critical current is 
^-independent, while at low temperatures it varies by 
a factor of two upon variation of <5$, see Fig. @] The os- 
cillations of the critical current as a function of the flux 
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through the normal region are reminiscent of the Fraun- 
hofer oscillations in a conventional Josephson junction, 24 
but the minima are not at zero and the periodicity is 
2(po rather than ipQ. These are characteristic signatures 
of a supercurrent carried by edge states rather than bulk 
states, 2 ^ 

The low-temperature supercurrent decays oc 1/L if 
the separation L of the NS interfaces is increased at 
constant width W. This is the characteristic decay of 
the spin-triplet proximity effect in a single transport 
channel— 



C. Meissner phase 

In the Meissner phase l m > A we may take an s- 
independent vector potential Aar. along each NS inter- 
face, given by Eq. (j2.10j) . If we also take s- independent 
parameters v c and «a, the two quantities A n and a n 
defined in Eqs. (|4.6I) and (14.71) are given by 



A„ = T A e z 



sin a r , 



, a n = TrWA AR /<p Q , (4.14) 



with ta = Wva I v 2 . (We kept the subscript n to allow for 
possibly different values of Aar at the two NS interfaces.) 

The zero-temperature limit (|4.11[) of the supercurrent 
then takes the form 



'o 

sinai sin 0:2 



(4.15) 



with the function T defined in Eq. (|4.13l) . 

The phase shift in the ^-dependence has disappeared, 
so now the supercurrent is an odd function of the flux $ 
through the ring, vanishing at $ = 0. Since dl/d$> > 
at $ = (for a.\ =0:2), the supercurrent is paramagnetic 
— in contrast with the usual diamagnetic Josephson ef- 
fect. Such a 7r-junction appears generically in the spin- 
triplet proximity effect. 22 The main effect of the phase a n 
accumulated by the vector potential along the NS inter- 
face is the reduction of the critical current by the factor 
(sinai sinai2)/(aia2) — the supercurrent vanishes if ot\ 
or ct2 is a (nonzero) integer multiple of tt. 

From Eq. (14. 15)) we conclude that the scaling of the 
critical current with the parameters of the Josephson 
junction is given in the Meissner phase by 



er 



A 



/2 
L rn 

w 2 



eu c (d/l so ) 2 (l m /C) s 



(4.16) 



with C = 2(L + W) the length of the perimeter of the 
normal region. (All coefficients of order unity are disre- 
garded in this scaling estimate, as well as any oscillatory 
dependence on W.) 



D. Mixed phase 

In the mixed phase £0 < l m < X the vector potential 
A n (s) along the NS interface depends on the configu- 
ration of vortices that have penetrated into the super- 
conductor. There is now a random phase shift of the 
supercurrent, both as a function of $ and The zero- 
temperature critical current reaches its maximal value 
7™ ax at <5<I> = — {ct\ + a2)<po/ 7T , given according to Eq. 
(EH} by 



27re, 
3^ 



(4.17) 



The 1/L 3 scaling with the separation of the NS interfaces 
is unchanged, but the scaling with the width W depends 
on the statistics of A n (s), which determines the statistics 
of A„ according to Eq. (|4.6j) . 

We have calculated the average of 7™ ax for a random 
variation of A n (s) as a function of s, with zero average 
and correlation length of order l m (the average separa- 
tion of vortices). The magnitude of the fluctuations is 
quantified by taking a piecewise constant A n (s) in each 
segment of length l m , drawn from a Gaussian distribu- 
tion with zero average and standard deviation axipo /Trl m 
with a of order unity. We have found that the average 
critical current in the mixed phase scales for W ~> l m as 



/7-max\ ^ eT A jm 
Vc I — a w 



d 2 ilw 



T 3 W 



c i2 r3 ' 

'so 

larger than in the Meissner phase by a factor W/l 



(4.18) 



V. COMPARISON WITH SPIN-SINGLET 
SUPERCURRENT 

A. Transfer matrix 

It is instructive to compare the results of the previous 
section for the spin-triplet supercurrent with the spin- 
singlet case considered by Ma and Zyuzini^ For that 
purpose we assume spin degeneracy in the 2D electron 
gas, neglecting Zeeman splitting or spin-orbit coupling. 
Electron-hole symmetry now relates excitations from op- 
posite spin bands, say an electron from the spin-up band 
and a hole from the spin-down band (or vice versa). 

The effective Hamiltonian of a spin-singlet edge chan- 
nel, to linear order in momentum, is 



H 



\{v c ,p~eA] A 

A* \{v c ,p + eA} 



(5.1) 



fully constrained by Hermiticity and the electron-hole 
symmetry requirement 



(JyH*a y 



H. (5.2) 
Choosing a gauge so that A is real we now have 

H = v c (p - eAa z ) + Aa x - \iv' c . (5.3) 
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The key difference with the effective Hamiltonian (12.5[) 
for a spin-triplet edge channel is that the coupling be- 
tween electrons and holes does not vanish at p = in the 
spin-singlet case. 

We now follow the same steps as in Sec. Ill El The 
particle current operator 



J = dHj dp = v c 

transforms H to 

H = J-^HJ- 1 / 2 =p- eAa z + (A/v c )a x , 
and produces the unitary transfer matrix 

e - Act 



(5.4) 



(5.5) 



M(s2, si) = V s exp 



S2 

I I ds 

»i 



eAa, 



(5.6) 

The transfer matrix no longer commutes with a z at e = 0, 
so there is no low-energy suppression of Andreev reflec- 
tion as in the spin-triplet case. The order parameter A 
equals Ao along the NS interface and zero along the nor- 
mal boundary. 



B. Meissner phase 

We consider the Meissner phase l m > A, with an s- 
independent vector potential A n along the interface with 
S„. Taking also an s-independcnt v c , we can evaluate Eq. 
(|5.6|) without the complications from operator ordering. 
The scattering matrix becomes 

S(e) = e ieT °e^' T /^ ^* + ' 5 */ 2 ^ cr ' I M2 

M n = exp[ieWA n cr z - i(A W/v c )a x ] , (5.8) 
with To = 2(L + W)/v c . The supercurrent follows from 
d x 

/(*) = ^ 2ksT ln det t 1 - s ( iuj p)] > ( 5J ) 

which differs from Eq. (|5.9p by a factor of 2 because of 
spin degeneracy of the edge channel in the spin-singlet 
case. 

Substitution of Eq. ([5~7) into Eq. (fO]) gives 



sin 2 /3 



AirkaT 

/($)= --^-sm{2^/ Vo ){W/^) 2 



x ^[cosh(wpro) +X] \ (5.10) 

p=0 



X = [cos(2^$/^o) - cos(7r<5$/^o)](W/e c ) : 
sin 213 

+ (tiWAar/Vo)— g— sm(ir6$/<p ) 

— cos2/3cos(7r<5$/<y9o)7 
(3 = ^(irWA AR /<p ) 2 + (W7£c) 2 . 



, sin 2 /3 



(5.11) 
(5.12) 



(For a compact expression, we took A\ = A^ = -Aar-) 
We defined the length £ c = Hv c /Aq, smaller than the 
superconducting coherence length £o = Hvf,s /Aq by a 
factor v c /vf,s- ln the point contact limit W — > con- 
sidered by Ma and Zyuzin our result (|5.10[) agrees with 
their finding (Eq. 13 of Ref. i). 

At zero temperature Eq. (|5.10j) evaluates to 



4p sin 2 ft 

/($)= sm(27r$/po)W£ c ) 2 ^ 

7TT 

1 / 1 — X 

x — . arctan 



VT^x^ 



VI- A 2 



(5.13) 



In contrast to the spin-singlet result (|4.11l) . the depen- 
dence of the supercurrent on the flux $ through the ring 
is strongly non-sinusoidal. The critical current oscillates 
both as a function of the flux S$ through the normal 
region and as a function of the width W of the NS inter- 
face. At high temperature only the oscillation with W 



sin 2 ft 

J($) = -8ek B Tsm(2w$/vo)(W/Z c ) 2b —^e-* kBTT °, 

(5.14) 

while the ^-dependence is now sinusoidal. 

Upon increasing the separation L of the NS inter- 
faces the spin-singlet supercurrent (|5.13[) in the low- 
temperature limit decays as 1/L. This is in contrast to 
the 1 /L 3 decay of the spin-triplet supercurrent (14.111) . In 
the high-temperature limit the supercurrent has the same 
exponential decay cx exp(— TrkgTTQ) in the spin-singlet 
and spin-triplet cases — only the pre-exponentials differ 
[cf. Eqs. (14~T0)) and ([5713]) ] . 

The spin-singlet supercurrent in the high-temperature 
regime has been studied also by Ishikawa and 
Fukuyama^i without taking the point contact limit 
W —> of Refs. 7.8J. We have not been able to reconcile 
their result with our Eq. (|5.14p . because only the length 
L of the normal boundaries enters into their exponential 
decay (rather than the sum L + W of the lengths of nor- 
mal and superconducting boundaries). The very recent 
study by Stone and Lin^ which also includes finite-W 
effects, still assumes W <C L so it does not distinguish 
between the two decay rates. 



C. Narrow-contact regime 



The full expression (|5.13l) for the zero-temperature 
spin-singlet supercurrent simplifies considerably in the 
narrow-contact regime W <C l m , when we may set 
Aar — »• 0. (This is the regime considered by Stone and 
LinJ£) Note that £ C /Z m ~ hcu c /A < 1, so W may still 
be large compared to £ c in the narrow-contact regime. 
As shown in Fig. [5j the current-phase relationship in the 
narrow-contact regime has a sawtooth-like shape, consis- 
tent with Ref. [H 

For reduced width w = W/£ c (modulo 7r) much less 
than unity the critical current exhibits resonant peaks (of 



8 



l - 



2 



-1 



-2 
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w/t 


r (mod 7r) = 

xS. — 01 
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0.5 



FIG. 5: Zero-temperature supercurrent as a function of the 
flux $ through the ring, for a flux 5$ through the normal 
region equal to an integer multiple of 2ipo — h/e. The green 
and blue curves (in units of (e/ro)| sinVF/^cl) are the spin- 
singlet result (|5.13[) for two values of W (in the narrow-contact 
regime W <C /,„, so with A ah — > 0). The black curve is the 
spin-triplet result (|4.1ip . plotted in units of ct\/tq. For the 
sake of comparison we also took the narrow-contact limit of 
the spin-triplet result, setting Qi,Q2 — > in Eq. (|4.1ip . 

W/i c (modTr) = - tt/2 Q.l 0.01 



C 1 A 




FIG. 6: Low-temperature critical current as a function of the 
flux 5$ through the normal region. The dashed curve is the 
spin-triplet result (|4.11|l . plotted in units of eT^/r§ (in the 
narrow-contact limit Qi,Q2 —¥ 0). The solid curves (in units 
of (e/ro)| sin iy/^c|) follow from the spin-singlet result (|5.13(l 
for three values of W in the narrow-contact regime. The res- 
onance at integer 5$/2(po peaks at I c = (27r/3)er^/ro in the 
spin-triplet case and at I c = (2e/ro)| sin VF/£ C | in the spin- 
singlet case. 



height 2ew/To) whenever 5$>/2ipo is an integer. (See Fig. 
El blue and red curves.) For w —> tt/2 the critical current 
I c = 2e/ro becomes (5<£>-independent (green line in Fig. 
[5]) — signifying the absence of Fraunhofer oscillations. 



two cases. Without spin-polarization we have the spin- 
singlet Joscphson effect considered earlier with low- 
temperature scaling 



srlet — ew c - 



(6.1) 



inversely proportional to the length C of the perimeter 
of the normal region. 

We have found that a spin-polarized Landau level can 
still carry a supercurrent. The low-temperature scaling 
of this spin-triplet Josephson effect is 



^c, triplet — ew c 



(W/l m )(d/l so ) 2 , 



(6.2) 



in the mixed phase with W 3> l m - 

For W ~ C the ratio of spin-triplet and spin-singlet 
critical currents is of order 



-^c, triplet /^c, singlet — (J>m/£) (d/lso*) 



(6.3) 



The spin-orbit scattering length in InAs is of order l so ~ 
100 nm, which could well be of the same order as the 
electrostatic length d (the smoothness of the potential 
step at the NS interface). The main reason for the rela- 
tive smallness of the spin-triplet supercurrent is then the 
factor l m /C. Since l m < 25 nm for B > IT, a submi- 
cron junction is needed for an observable effect. As we 
have shown, the spin-triplet Josephson effect has unusual 
features, including a paramagnetic, rather than diamag- 
netic, current-phase relationship, and Fraunhofer oscilla- 
tions which have a h/e rather than h/2e periodicity. 

For the purpose of comparison with the spin-singlet 
Josephson effect, we have performed an analysis that 
goes beyond earlier work on that problemj£&i£ in par- 
ticular with regard to the Fraunhofer oscillations. We 
have found a remarkable dependence of the amplitude of 
the Fraunhofer oscillations on the relative magnitude of 
the junction width W and an effective coherence length 
£ c . For W/^ c = tt/2 (mod it) the Fraunhofer oscillations 
vanish alltogether, see Fig. [51 

These spin-singlet results may well be of relevance also 
for graphene, which is an attractive alternative to InAs 
in the search for the coexistence of the Josephson and 
quantum Hall effects. The results obtained here would 
apply if W is larger than the intervalley scattering length. 
For smaller W the valley-selectivity of the edge states 
enters, along the lines described in Ref. HH. 



VI. CONCLUSION 

In conclusion, we have analysed the Josephson ef- 
fect in the lowest Landau level, both with and without 
spin-polarization. The critical current scales differently 
with the parameters of the Josephson junction in these 
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Appendix A: Andreev-Rashba edge states 

The theory of Andreev edge states, produced by the in- 
terplay of cyclotron motion and Andreev reflection, has 
been developed by Ziilicke and collaborators ! 4 ! 20 ' 28 Here 
we include the interplay with Rashba spin-orbit interac- 
tion, in the spin-polarized regime where Andreev reflec- 
tion can only occur because of the Rashba effect. 

The theory is complicated by the fact that we are deep 
in the quantum mechanical regime, with only one occu- 
pied Landau level, and cannot make the semiclassical ap- 
proximation of large Landau level index made in earlier 
work i 4 ' 20 ' 28 ' 29 Since the Fermi energy in the normal metal 
is small compared to the superconducting gap, we can 
also not make the usual Andreev approximation (match- 
ing wave amplitudes without matching derivatives). We 
keep the theory tractable analytically by treating the 
spin-orbit interaction perturbatively. 

The goal of our analysis of the Andreev-Rashba edge 
states is to arrive at a microscopic derivation of the pa- 
rameters that enter into the effective edge state Hamilto- 
nian (|2 .3[) , on which our theory of the spin-triplet Joseph- 
son effect is based. 



1. Bogoliubov-De Gennes equation 

We start from the Bogoliubov-De Gennes (BdG) equa- 
tion 

fH -E F Ar y\(f)=e(f), (Al) 

for quasiparticle excitations consisting of an electron 
spinor ip e = (u + ,u_) and a hole spinor iph = (t; + ,w_). 
The label ± indicates the spin band and the Pauli matrix 
T y acts on the spin degree of freedom. The pair potential 
A of a spin-singlet superconductor couples electron and 
hole excitations in opposite spin bands. Electron-hole 
symmetry is expressed by a x H*a x = —H. 
The single-particle Hamiltonian 

H = (2m)- 1 (p - eAf + V + \g m B • r + H R , (A2) 
Hr = a(p y - eAy)T x - a(p x - eA x )r y , (A3) 

contains the kinetic energy, potential energy, Zeeman en- 
ergy, and Rashba spin-orbit interaction. We consider a 



translationally invariant NS interface at y = 0, with vec- 
tor potential A = A(y)x, magnetic field B — —A'(y)z, 
electrostatic potential V — V(y), and pair potential 
A = A(y). The effective mass m, effective gyromagnetic 
factor g, and Rashba coefficient a arc taken spatially uni- 
form (otherwise also derivatives of m and a would have 
to enter in the Hamiltonian, to preserve Hermiticity). 

Parallel momentum p x = p is conserved, for states 
oc e lpx . The y-dependence of the wave functions is deter- 
mined by 

- \g^A'{y)T z + H R , (A4) 
H R = -iaT x -^--a[p-eA(y)}Ty. (A5) 

In this basis the operators Hq,Hq in the BdG Hamilto- 
nian should be replaced by H (p), H^—p). 

The NS interface is at y = 0, with the superconduc- 
tor in the region y < 0. In the simplest model for the 
interface we take a step function both for the pair po- 
tential, A(y) = Ao#(— y), and for the electrostatic po- 
tential, V(y) = -V 0(-y) with V > 0. (The function 
9{y) equals 1 for y > and for y < 0.) Smoothing of 
the interface is important, and will be considered at the 
end of the Appendix. We assume that we are deep in 
the Meissner phase, l m ^ A, so that we may neglect the 
penetration of the magnetic field in the superconductor. 
In the gauge where Ao is real, the vector potential is then 
given simply by A(y) = —yB0(y). 

We will first solve the eigenvalue problem to zeroth 
order for H R = 0, and then include the Rashba spin- 
orbit interaction to lowest order as a perturbation. 



2. Solution without the Rashba effect 

a. Eigenstates in S 

In S (for y < 0) the BdG Hamiltonian with H R = is 
given by 



/Kp) 






*A 



Kd v 





H(p) - nd\ 
^A 






zA 

-H{p) + K&1 





-zA 



-fj,(p) + nd 2 y) 



(A6) 



with fi(p) = p 2 /2m — Vq — Ep, k — (2m) 1 , and d y = There are four eigenstates Xs,±w±(y) of H$ for < 

d/dy. e < Aq (decaying for y — s- -co), with 



w ± (y) = e^y, 7 ±^±M/A 2 -e 2 , (A7) 
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nq 2 ± = -(i(p) ± i^A 2 -e 2 , Im q ± < 0, (A8) 

(A9) 







( ° 


V 








, Xi,± = 




iA J 




V o 



Xt,± 



For A <C Ep + Vq = Ep,s = P% s /2m we may approxi- 
mate 



«± = Tq(p) - ^y\/A^-£ 2 , g(p) = ^fp FiS -p 2 - 

(A10) 



6. Eigenstates in N 



In N (for y > 0) we have, again for TJr = 0, 



H 



N 



(U{p,y)~Kd 2 y +n+ 

U{p,y)- K d 2 y +H- 

-U(-p,y) + Ktf- fi+ 











with /z± = ±±g/j, B B - E F and U(p,y) = (p + eBy) 2 /2 



rn . 







-U(-p, y) + Kd 2 .- n_J 



(All) 



The differential equation 

[U (p, y) - nd 2 y + fj,±]<f>(y) = e</>(y), 



(A12) 



with <j)(y) — > for y — > — oo is solved by a parabolic 
cylinder function 



V2[^-+pl v 

'm, 



. (A13) 



The normalization constant C^ p — 0(lm^ 2 ) is deter- 
mined by 

poo 

/ 2 ± (e,p,y)dy = l. (A14) 
Jo 



The parabolic cylinder function lA{—v^y) has no nodes 
as a function of y for v < 1/2 and only a single node for 
1/2 < i/ < 3/2. 

The four eigenstates of ifjv are constructed in terms 
of the functions <j>± , 



V'et = 



/0+(e,p,»)\ / 



-(-e,-p,y)J 

(A15) 



c. Matching at the NS interface 

We construct two independent superpositions of basis 
states in N and S, 

= %) [ai^, t (y) + 6i^e,4.(l/)] 

+ 0{-v) [ciXU+e lq+v + d lX u-e lq - y ] , (A16a) 

*2(y) = %) [a 2 ^ e ,t(y) + 62^4 (y)] 

+ [c 2 XT,+e <9+!/ + rf 2 Xt,-e i9 - !/ ] • (A16b) 

We choose e such that 

(s - n+)/hLJ c < \ < (s - H-)/tvbJ c < f. (A17) 

In this range the equation (j) + (s,p, 0) = has no solution 
while the equation <j>_(e,p, 0) = has a single solution 
p = T>(e). As we will see, this is the branch of the dis- 
persion relation with wave function while another 
branch, with wave function ^2, is given by p = —V(—e). 

Continuity of and d^i/dy at y = gives four equa- 
tions for the coefficients oti, b\, c\, d±, 

a\(f> + {-e, -p,0) = ci7_ + dij+, (A18a) 

b 1 <t>-(e,p,0)=i&o(ci+d 1 ), (A18b) 

Oi0' + (-e, -p, 0) = 2<?+ci7_ + 2(j_rfi7 + , (A18c) 

6i0'_(e,p,O) = -A (<z + c 1 + g _d 1 ), (A18d) 

with </>± = d(j)±/dy. The solution satisfies 

ci _ 7+ q-(t>+(-e,-p,0) + i<f>' + {-e,-p,0) 



di 7- q+(p + (-e,-p,0) + i<p' + (-e, -p,0) 



7+g- 
7-9+ 



[1 + 0(A F , s /U] , 



(A19) 
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since <f>' + is smaller than q±<fi+ by a factor XF,s/lm "C 1 
(with Af.s = 27r/fc F .5). [Here we have used that <j> + does 
not vanish for e in the range (|A17[) .] 

Similarly, for ^ 2 we have the matching conditions 

a 2 <f>+(e 7 p,0) = c 2 7+ + d 2 7_, (A20a) 

& 2 0_(-£,-p,O)-iA o (c 2 +d 2 ), (A20b) 

a 2 4>' + (e,p, 0) = i<7+c 2 7 + + iq-d 2 j-, (A20c) 

6a^_(-e,-p,0) = -A (g+c 2 + g_da), (A20d) 

with solution 

c 2 7_ g_</> + (e,p,0) + i(p' + (e,p,0) 



d 2 -y+ q + <j) + (s,p,0) +i<p' + (e,p,0) 



i-q- 

7+94 



[1 + 0{X F , s /l m )] 



(A21) 



The normalization requirement gives one more equa- 
tion for each set of coefficients, 



K| 2 + |6„| 2 + 



These results provide the numerical coefficients for v c and 
A AR in Eqs. (JSU) and (|2T0]t . 

Concerning the coefficient c m , we note that, as required 
by Eq. (|A17[) , the value of v_ at the Fermi level is in the 
range 1/2 < v_ < 3/2. The ratio g^B/hio c — gm/2mo 
(with mo the free electron mass) is typically much smaller 
than unity, so v- « 1/2 will be close to the lower end of 
this range and c m s» 0.88. The dispersion curves in Fig. 
[2] are plotted for V- = 1/2. 



e. Eigenstates 

From the matching conditions we determine the coef- 
ficients of the zeroth order eigenstates, 



a± 2iAo 

di <fi + (—Ep,—p,Q) 



a 2 _ -2iA _ 

d 2 (j> + (-£- p ,p,0) 

(A30a) 



K| 2 ) = l. (A22) di #_(e p> p,0) 



b x _ ~2A q(p) b 2 _ 2A q{p) _ 

= A 1' T ~ ITT- — = A 2- 



d 2 4>'_(e- p ,-p,0) 



(A30b) 



It follows that 



Since 



7+9- 



d. Dispersion relation 



l + 0( £ /A o ) + 0(A F ,s/£o), (A23) 



7-9+ 

we may approximate 

-r = -1 = -r for £ < A and A F s < £o,Jm- (A24) 
ai a 2 

Eqs. (|Al8b|) and (|A20bD then give the dispersion relations 



Ei(p) = 2? in » = £ p , for (A25a) 
E2 (p) = -P inv (-p) = -£- P , for * 2 , (A25b) 

with e p determined by the equation 

<P-{e p ,p,0) = i)^U ^~ Ep ,V2pl m = 0. (A26) 

The dispersion relation of the two modes is plotted in 
Fig. [2] For small p it is approximately linear, given by 



£ P = v c (p- eA AR ), 



(A27) 



with the definitions 



v c = 1.14 l m oj c , ew c AAR = {v- - |)w c , 



• A AR = 0.88 - \)l m B = -c m l m B, (A28) 

E ? + y^ B . (A29) 
nio c 



a n /b n = 0(X F ,s/l m ) < 1. 



(A31) 



This means that VPi and V E' 2 in the normal region have 
most of their weight in the spin-down band, so is 
predominantly an electron state and ^ 2 is predominantly 
a hole state. 

The normalization condition (|A22I) simplifies to 

|6„| 2 + ^>K| 2 = 1, Y = 2q{p)A /m. (A32) 



Together with Eq. (|A30[) this determines all coefficients 
(up to an overall phase factor), 

a n = Y n (X 2 + Yq)- 1 ' 2 , b n = X n (X 2 + Y )- 1/2 , 
d n = -c n = (Xl + Y )- 1/2 . (A33) 

Because 4>'_(s p ,p,0) = 0(lm 3 ^ 2 ), we can estimate 



— = O 

X 2 



/ 3 



« 1, 



(A34) 



since we work in the regime where l m is large both com- 
pared to Xf,s an d compared to fo- We may therefore 
neglect Y relative to A 2 . 



3. Inclusion of the Rashba effect 



We include the Rashba Hamiltonian 



SH = 



Hr, 





-Hi 



(A35) 
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as a perturbation of the BdG Hamiltonian. To lowest 
order in this perturbation we need the matrix elements 
of SH in the basis of unperturbed eigenstates ^i,^> 2 . 
Since (fy n | SH \ $> n ) = 0, there is only a single matrix 
element (* 2 | SH | *i) = (*i | SH | <f 2 )* to consider. 
We calculate separately the contributions to this matrix 
element from the superconducting and normal regions. 



a. Matrix element in S 

The Rashba Hamiltonian in the superconducting re- 
gion is 



SH S = 
Note that 



_ (—ioLT x d v — apr y 

' —iaT x d y + apr y 



(Xt,± I SH S | Xl,±)s = 0, 



(A36) 



(A37) 



where (■ ■ - )s indicates integration over the superconduct- 
ing region y < 0. The matrix element becomes 



(# 2 | SH S | = iaA y - s 2 
x [<^d x {l + ip/qj) - Cl d* 2 (l + ip/q+)} 

With the help of the approximation 
c 2 d\ 



c\d% 



(A38) 



(A39) 



this gives, for e <C A , 

(* a | SH S | *i)s = aA 2 oCl d* 2 (p/q + - p/q_). (A40) 

Since q± = TPF,s{l + 0(\ F ,s/£o)+0(p/p F ) 2 }, we may 
further approximate 



(* 2 | SH S | = -2aA 2 oCl dl 



V = 2aA 2 Q p 
Pf,s Pf,s\XiX2\ 



ap 



(kF,slm) 3 



*s(p), 



(A41) 



where we have used Eqs. (|A33[) and (|A34|) . and we have 
introduced a dimensionless even function of p, 



*s(p) = |4|^_(e-p ) -p ) O)0 , _(ep,p ) O)|. 



(A42) 



See Fig.[7]for a plot of p<&s{p), which is an approximately 
linear function of p, given for small p by 

Up*s(p) = 1-13 ImP + 0{l mP ) 2 . (A43) 



b. Matrix element in N 



The Rashba Hamiltonian in the normal region is 



The matrix element is 



SH N 



_ f—iar x d y — a(eBy + p)r y 

-iar x d y - a(eBy - p)r y 



(A44) 



(* 2 | SH N | *i)jv = iaalbi{<l> + (-E-p,p,y) \ -d y + eBy+p \ cj)^(e p ,p,y)) N 

+ iablax{(t>-{e- p ,-p,y) \ -d y - eBy+p | (f) + (-e p ,-p,y)) N , 



(A45) 



with the coefficients given by Eq. (|A33|1 . 



in terms of a dimensionless even function of p, 



a 2°l = - , v v , , Ml 



XxX 2 



\x 1 x 2 \ 



(A46) 



The matrix element can be written in the form 



(* 2 I SH N | *i) 



N 



ap 



F,S l m 



-*at(p), 



(A47) 



$jv(p) = *(p)+*(-p), 

1 W_{e_ p ,-p,Q)4>'_{e p ,pM 



(A48a) 



P <l>'-(£p,P,0)<l>+(-£-p,P,0) 



x (<f>+(-£-p,P,y) I -j-+plm~ l m d y I 4>-{e p ,p,y)) N . 

''771 

(A48b) 

In Fig. [7] we have also plotted pQ n (p) ■ The p-dependence 
is approximately linear, given for small p by 



lmP$N{p) = C N l rn p + 0(l m p) 2 



(A49) 



13 



e< l 

*«| 

~"-l 

-2 



1 1 










' N 






* --f 






1 



-1 1 

P [fl/l m } 

FIG. 7: Plot of the functions p$s(p) and p$jv(p), which 
determine the contribution to the Rashba matrix elements 
(|A41[I and (IA47[) from the superconducting and normal 
region, respectively. The curves are calculated from Eqs. 
(|A42|) and (|A48|) . for Fermi energy E F = ifoj c . These 
two functions are of the same order of magnitude, but the 
contribution to the Rashba matrix element from S has an 
additional prefactor (kF,sl m )~ 2 , so it is much smaller than 
the contribution from N. 



The coefficient cat is a function of E F /hio c , of order unity. 
For E F — ^huj c (Fermi level half-way between the splin- 
split lowest Landau level) one has c/v = 1.98. 



Andreev-Rashba edge states at an abrupt NS 
interface 



From Eqs. (|A4l"j) and (lA47l) we find that the matrix 
element of the Rashba Hamiltonian in the unperturbed 
basis is 



(* 2 I SH I *0 



ap 



[<S> N (p) + (kp S l m r 2 <S>s(p)]- 



kF,slr, 

(A50) 

Since both functions $jv(p) and <&s(p) are of order unity 
for p of order l/l m , the effect of spin-orbit coupling in 
the superconductor on the Andreev-Rashba edge states 
is weaker by a factor l/(kF,slm) 2 — f^c/Ep,s than the 
effect of spin-orbit coupling in the normal region. We 
therefore arrive at the final result for the Rashba matrix 
element, 



(* 2 I SH I *i) 



ap 



-$jv(p)- 



(A51) 



To first order in the Rashba coefficient a, the BdG 
Hamiltonian in the unperturbed basis is a 2 x 2 matrix 
T-L with elements 



H = 



(ap/k F ,slm)$N(p) 



(ap/k F , s l m )$N(p) 



(A52) 



The matrix elements have an approximately linear p- 
depcndcncc, 



v c (p-eA A R) v A p 

vaP v c (p + eA AR ) 



(A53) 



with coefficients v c and Aar given by Eq. (|A28|) . The 
coefficient va follows from Eq. (|A49|) . 



v A = c N 



kF,slm 



'tF.sho 



(A54) 



in terms of the spin-orbit scattering length l so = h 2 /ma. 

The dispersion relation of the Andreev-Rashba edge 
states, to second order in the Rashba coefficient a, is 
given by 



e± = v c p ± yj (ev c A AR ) 2 + {vap) 2 , 



(A55) 



with the + sign for the electron-like mode ^Pi and the 
sign for the hole-like mode ^2- 



Andreev-Rashba edge states at a smooth NS 
interface 



So far we have taken an abrupt model for the NS in- 
terface, with a step function both in the pair potential 
(from to Ao) and in the electrostatic potential (from 
to — Vo). We now turn to the more realistic model 
of a smooth interface. Since Ao <C Ef,s we do not ex- 
pect the abruptness of the pair potential step to have 
significant consequences, so we keep the step function 
A(y) = Ao6(-y). 

The situation is different for the electrostatic potential 
step, which enforces normal reflections at the expense 
of Andreev reflections. We therefore broaden the step in 
V(y) over a distance d, such that V(y) = — Vo for y < —d 
and V(y) — for y > 0. The abrupt limit corresponds 
to d ~ 1/kp.s- We now take d larger, but still small 
compared to £o- 



a. Eigenstates m S 

The eigenstates in N are unaffected by the smoothing 
for y < 0. The eigenstates in S are given by Xs,±w±(y), 
with the same spinor Xs.± defined in Eq. (|A9I) and a 
spatial profile w±(y) determined by 



d 2 y w ± (y) 



2m 



-Kp,y)± i \f£$ 



l*(p,y) =p 2 /2m + V(y) -E 



F- 



w±(y), 
(A56) 
(A57) 



Since we assume d <C £ o w e may solve the scatter- 
ing by the potential step independently of the reflection 
from the pair potential. The wave vector (in the limit 
e — > 0) changes from k(p) 
y/2m(E F 



k'(p) = y]2m{E F + Vo 
lutions a + e lky + a-e~ lky at y 
wave solutions a' + e lk v 
scattering matrix, 



\j2mEp — P 2 at y = to 
p 2 at y — — d. Plane wave so- 
are related to plane 
at y = — d by a unitary 



r t 
t' r 1 



(A58) 



14 



The solution w±(y) corresponds to setting aL — 0. We 
thus obtain 



1 tu'_(0) 



(0) 



r + 1 iw_(0) 



1 



(A59) 



with w+ 



dw±/dy. The complex conjugation appears 
as a result of inversion of the scattering matrix, but it 
can be ignored because the reflection amplitude r is real 
for kd < 1. 



b. Matching at the NS interface 

Matching of the eigenstates in N to those in S at 
y = proceeds entirely as in the case of the abrupt in- 
terface, with q± replaced by the logarithmic derivative 
w±(0)/zw±(0). The result (1A"30|) for the matching coef- 
ficients changes simply by the replacement of q(p) by 



IcsiP) = — 7^ = k{P)- 



iw-(0) 

The Fermi wave vector kp.s = <?(0) is replaced by 

1 -r 



fccff = geff(O) = k(0)- 



+ r 



(A60) 



(A61) 



The reflection amplitude r is related by r = — yl — T 
to the over-barrier transmission probability T. Since 
d <C Co and fc ~ ^ 1/Coj we necessarily have 

fee? ~ (i/irn ^ 1- One may then expand 



T = Cbarriorfcrf + 0{kd)' 



(A62) 



with Cbarricr a numerical coefficient of order unity. The 
wave vector k e g takes the form 



4fc(0) 



Cba 



(A63) 



The value of Cbarrier depends on the shape of the bar- 
rier. As an example, we take the Woods-Saxon step 



V(y) = -V [l + e {v+Vo)/d ] \ 



(A64) 



with yo ^> d (so that the potential is essentially zero for 
y > 0). The transmission probability is^ 



T = 1 



sinh 2 [7rd(fc- k 1 )] 



sinh 2 [7T(i(fc + fc')] 
= 47rfcdcotanh(7rfc'd) + O(kd) 2 . 



(A65) 



For a step which is smooth on the scale of Xp,s ( so k'd ~ 
d/XF,s 3> 1) we arrive at Eq. (j A62|) with Cbarricr = 47r, 
hence fc e ff = l/nd. In the opposite regime k'd — >■ of an 
abrupt potential step we have T = 4fc/fc', hence k e g = 
fc'(0) — as it should be. 

c. Dispersion relation 

The zeroth order dispersion relation, which is indepen- 
dent of q(j>), remains unchanged, still given by Eqs. (|A25|) 
and (|A26[) , and also the velocity v c remains given by Eq. 
(|X28]l . 



The expression (|A51j) for the Rashba matrix element 
in the unperturbed basis is changed into 



(* 2 I SH I *i) = 



ap 
k c gl„ 



■*at(p), 



(A66) 



with the same function $jv a s for the abrupt interface 
(see Fig. [7]). Once again, the dominant contribution to 
the matrix element comes from the normal region, with 
the contribution from the superconducting region smaller 
by a factor (k c sl m y 2 — (d/l m ) 2 . 

The BdG Hamiltonian in the unperturbed basis still 
has the form (|A53[) . the only difference appearing in the 
coefficient ua ■ Instead of Eq. (|A54|) for the abrupt inter- 
face it is now given by 



va = c N 



Vgd 



(A67) 



This is the result (|2.8p used in the analysis of the spin- 
triplet Josephson effect. 
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